Nonextensive Thermostatistics and the //-Theorem 
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The kinetic foundations of Tsallis' nonextensive thermostatistics are investigated through Boltz- 
mann's transport equation approach. Our analysis follows from a nonextensive generalization of the 
"molecular chaos hypothesis". For q > 0, the g-transport equation satisfies an //-theorem based on 
Tsallis entropy. It is also proved that the collisional equilibrium is given by Tsallis' g-nonextensive 
velocity distribution. 
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In 1988 Tsallis proposed a striking generalization of 
the Boltzmann-Gibbs entropy functional given by [l], 



qPi 



(1) 



where k is Boltzmann's constant, pi is the probability of 
the i-th microstate, the parameter q is any real number, 
and the y-logarithmic function is defined as 



hx q f = (l-q)- L (f 



-l [f l-q 
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For future reference it is convenient to introduce the q- 
cxponential function e q (f), which is defined by 



e q (f) = [1 + (1 - q)f] 
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(3) 



if 1 + (1 — q)f > and by e q (f) — otherwise. These 
q-functions satisfy, for f,g>0, the identities (see || for 
a thorough discussion) 



e g( ln <?/)=/ i 
ln g / + ln g g = \n q fg 



(q-l)(\n q f)(\n q g) 



(4) 



When q — > 1 all the above expressions reproduce those 
verified by the usual elementary functions and Tsallis' en- 
tropy reduces to the standard logarithmic one, namely: 
Si = —k'^2 i pi\iipi. The most distinctive trait of S q is 
its pseudoadditivity. Given two systems A and B in- 
dependent in the sense of factorizability of the (joint) 
microstate probabilities, the Tsallis entropy of the com- 
posite system A®B verifies S q {A®B) = S q (A) + S q (B) + 
(1 — q)S q (A)S q (B). Hence, |1 — q\ quantifies the lack of 
extensivity of S q . 

The q-thermostatistics associated with S q [pl^j is nowa- 
days being hailed as the possible basis of a theoreti- 
cal framework appropriate to dealwith nonextensive set- 
tings. There is a growing body of evidence suggesting 
that S q provides a convenient frame for the thermosta- 



fHlllJlflL, sucn as tne ve l oc ity distribution of galaxy 
clusters || , Landau damping in plasmas 0, superdif- 
fusion phenomena || and, more generally, systems ex- 
hibiting a nonextensive thermodynamic behaviour due, 
for instance, to long range interactions [0J|]. It is worth 
to stress that some of the aforementioned developments 
involve a quantitative agreement between experimental 
data and theoretical models based on Tsallis' thermo- 
statistics j§|^||Q]. For instance, it was experimentally 
found that pure electron plasmas in Penning traps re- 
lax to metaestable states whose radial density profiles 
do not maximize Boltzmann-Gibbs entropy. However, 
Boghosian showed that the observed profiles are well de- 
scribed by Tsallis' thermostatistics with q close to 1/2 
H . Beck's recent treatment of fully developed turbulent 
flows constitutes another interesting example [pj. Based 
on Tsallis formalism with q ^ 1, Beck calculated the 
probability density functions of velocity differences de- 
pending on distance and Reynolds number, as well as the 
concomitant scalling exponents, finding good agreement 
with turbulence experiments. 

A large portion of the experimental evidence support- 
ing Tsallis proposal involves the non-Maxwellian (power- 
law) velocity q-distribution associated with Tsallis' gen- 
eralized canonical ensemble approach to the classical N- 
body problem. This equilibrium g-distribution may be 
derived through a simple nonextensive generalization of 
the Maxwell ansatz [[)) . Alternatively, it can be obtained 
maximizing Tsallis entropy under the constraints im- 
posed by normalization and the energy mean value p0[ , 
a procedure closely related to Jaynes information theory 
formulation of statistical mechanics |jll|,^2|. So far, most 
theoretical studies on Tsallis thermostatistics have been 
developed on the basis of the maximum entropy principle 
p], [13 14 ll|. As widely known, this approach to statis- 
tical ensembles was, in the case of standard statistical 
mechanics, historically the latest one to appear. Even 



tistical analysis of many physical systems and processes when Gibbs ' introduced his ensemble approach, the ki- 
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netic foundations of statistical mechanics were already 
well developed. On the light of this, it is not unreason- 
able to expect that a systematic exploration of the kinetic 
aspects of Tsallis' thermostatistics may be crucial to il- 
luminate its foundations as well as to achieve a better 
understanding of its physical applications. 

In this framework, the aim of this letter is to obtain 
the equilibrium velocity g-distribution from a slight gen- 
eralization of the kinetic Boltzmann ii-theorem. The 
whole argument follows simply by modifying the molec- 
ular chaos hypothesis, as originally advanced by Boltz- 
mann, and generalizing the local entropy expression in 
according to Tsallis proposal. 

The statistical content of Boltzmann's kinetic theory 
relies on two main ingredients Jl6|,[l7]]. The first one is 
a specific functional form for the local entropy, which is 
expressed by Boltmann's logarithmic measure 



H[f] 



k / f(x, v, t) In fix, v, t) d 3 v 



(5) 



The second one is the celebrated hypothesis of molec- 
ular chaos ("Stosszahlansatz"), which is tantamount to 
assume the factorizability of the joint distribution asso- 
ciated with two colliding molecules 



f{Xl,Vl,X2, ifa, t) = /(xi, Vx,t)f (x 2 ,v 2 ,t) 



(6) 



These two statistical assumptions are inextricably inter- 
twined. Therefore, if one adopts a generalized nonex- 
tensive entropic measure, a consistent generalization of 
the "Stosszahlansatz" hypothesis should also be imple- 
mented. Different choices for the collision term in the 
kinetic equation (which, in turn, is determined by the 
"Stosszahlansatz" ) lead to different forms for the entropic 
functional exhibiting a time derivative with definite sign. 
As a consequence, the form of the entropic functional 
behaving monotonically with time (and consequently ad- 
mitting an _ff-theorem) depends upon the form of the 
collisional term appearing in the kinetic equation. 

Historically, the basic assumption (|^) (sometimes re- 
ferred to as "Maxwell's ansatz") has generated a lot of 
controversy. The fundamental role played by this hy- 
pothesis was first realized by S.H. Burbury E^] in 1894. 
The precise characterization of the conditions of its ap- 
plicability has been an important conceptual problem of 
theoretical physics ever since [[l9| . The physical meaning 
of equation (|J) is that colliding molecules are uncorre- 
cted. The irreversible behaviour of Boltzmann's equa- 
tion can be traced back to this assumption. It is clearly a 
time asymmetric hypothesis, since molecules assumed to 
be uncorrelated before a collision certainly become corre- 
lated after the collision has taken place pfl] . 

Although very plausible, and endowed with an intu- 
itively clear statistical meaning, Boltzmann's particular 
expression (^) for the molecular chaos hypothesis can not 
be deduced from first principles. By no means it is an in- 
escapable consequence of classical mechanics. Boltzmann 



himself accepted that the hypothesis of molecular chaos 
was needed in order to obtain irreversibility. Further, 
He also admitted that the hypothesis may not always be 
valid for real gases, especially at high densities [ p^|pl] ] . 

In what follows we introduce a consistent general- 
ization of this hypothesis in accordance with Tsallis' 
nonextensive formalism. We remark that equation (j^) 
implies that the logarithm of the joint distribution 
f(xi, Vi,X2, V2,t) is equal to the sum of two terms, each 
one involving only the marginal distribution associated 
with one of the colliding molecules. Our generalized hy- 
pothesis is to assume that a power of the joint distribu- 
tion [instead of the logarithm) is equal to the sum of two 
terms, each one depending on just one of the colliding 
molecules. By recourse to the g-generalization of the log- 
arithm function, this condition can be formulated in a 
way that recovers the standard hypothesis of molecular 
chaos as a limit case. 

Let us now consider a spatially homogeneous gas of 
N hard-sphere particles of mass m and diameter s, un- 
der the action of an external force F, and enclosed in 
a volume V. The state of a non-relativistic gas is ki- 
netically characterized by the one-particle distribution 
function f(x,v,t). The quantity f{x,v,t)d 3 xd 3 v gives, 
at each time t, the number of particles in the volume 
element d 3 xd 3 v around the particle position x and veloc- 
ity v. In principle, this distribution function verifies the 
q-nonextensive Boltzmann equation 



+ v ' ^ H F ' = C i (/) 

at ox m ov 



(7) 



where C„ denotes the q-collisional term. The left-hand- 
side of Mj is just the total time derivative of the distri- 
bution function. Hence, nonextensivity effects can be in- 
corporated only through the collisional term. Naturally, 
C q (f) may be calculated in accordance with the laws of 
elastic collisions. Its specific structure must lead to the 
standard result in the limit q — ► 1. We also make the 
following assumptions: (i) Only binary collisions occur 
in the gas; (ii) C q (f) is a local function of the slow vary- 
ing distribution function; (iii) C q (f) is consistent with 
the energy, momentum and particle number conservation 
laws. 

Our main goal is now to show that the generalized 
collisional term C q {f) leads to a nonnegative expression 
for the time derivative of the g-entropy, and that it does 
not vanish unless the distribution function assumes the 
equilibrium form associated with q-Maxwellian gas ||. 
Now, following standard lines we define 



C q {f) = V / \V ■ S\R q dtod 3 v 1 



(8) 



where d 3 v\ stands for the volume element in velocity 
space, V = v% — v is the relative velocity before colli- 
sion, e denotes an arbitrary unit vector, and duj is an 
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elementary solid angle such that s 2 duj is the area of the 
"collision cylinder" (for details on the collision's geome- 
try see Refs. The quantity R q (f,f) is a differ- 
ence of two correlation functions (just before and after 
collision), which are assumed to satisfy a g-generalized 
form of the molecular chaos hypothesis. In the present 
q-nonextensive scenario we will assume that 



R g {f, f) = e q {f q ~ X hx q f + fT 1 In, /() 



,(/«-Mn 9 / + /r 1 ln 9 / 1 ) 



where primes refer to the distribution function after col- 
lision. When q — ► 1 equation reduces to Boltzmann's 
molecular chaos hypothesis 

lim R q = R = ff\ - fh . (10) 
For the local entropy we adopt Tsallis expression, 



H q = -kJ f q \n q fd 3 v 



(11) 



which reduces to the standard Boltzmann measure (||) 
for q = 1. Now, we first take the partial time derivative 
of the above expression 



dH q 

at 



-k J[ q p-Hu q f + l]^d 3 v . (12) 



As one may check, by inserting the generalized Boltz- 
mann equation (Q) into (|l^), and using (||), expression 
( |l2| ) can be rewritten as a balance equation 



dH, 



dt 



^ + V -S q = G q (f,t) , 



where the g-entropy flux vector S q associated with H q is 
defined by 



S q = -k I vf q In, fd 3 v 



(13) 



and the source term G q reads 



G q = -^- J \V-e](l + qf q - 1 \n q f)R q dcjd 3 v 1 d 3 v . 

In order to rewrite G q in a more symmetrical form some 
elementary operations must be done in the above expres- 
sion. Following standard lines [lq] , we first notice that 
interchanging v and v{ does not affect the value of the 
integral. This happens because the magnitude of the rel- 
ative velocity vector and the scattering cross section are 
invariants. Similarly, the value of G q is not altered if we 
integrate with respect to the variables v' and v[ (we recall 
that d 3 vxd 3 v = d 3 v[d 3 v'). Note that this step requires 
the change of sign of R q (inverse collision). Implementing 
these operations and symmetrizing the resulting expres- 
sion, one may show that the source term can be written 
as 



G q (f, t) = -^J\V- eKqfr 1 ln q fx + 3/ 9 " 1 / 

-qfT 1 K fx - qf' 9 ' 1 In, f)R q dwd 3 Vxd 3 v . (14) 

Making now the transformation f q hx q f = hx q * f, 
where q* = 2 — q, and rearranging terms we find 



G q — 



ks 2 q 



j\V- e1(hv /' + h V f\ - ln r / - h V fx) 



(9) [e q (ln r f + In,, f\) - e q (ln q * f + ln g , fx)}dud 3 vxd 3 " 



Note that the integrand in the above equation is never 
negative, because (ln g * /' + ln g » f[ — ln g * / — ln g * /i) and 
[e g (ln g . /' + In,. f\) - e q (ln q * f + ln g « fx)] always have 
the same signs. Therefore, for positive values of q, we 



obtain the 77 g -theorem 



dt 



V-S q = G q (r,t)>0. 



(15) 



This inequality states that the g-entropy source must be 
positive or zero, thereby furnishing a kinetic argument 
for the second law of thermodynamics in the framework 
of Tsallis' nonextensive formalism. However, our argu- 
ment does not constitute a kinetic proof of the second 
law. As happens with the standard Boltzmann equa- 
tion, our generalization can not be obtained only from 
the Hamiltonian equations of motion. Specific statistical 
assumptions are also needed. 

When q < the entropy of a given volume element 
decreases with time. Consequently, it seems that within 
the present context, and according to the second law of 
thermodynamics, the parameter q should be restricted to 
positive values |22|. Notice also that the entropy does 
not change with time if q = 0. Similar results were 
previously obtained using the master equation and the 
relaxation time approximation p3| ]. Naturally, Tsallis' 
g-parameter may be further restricted by other physical 
requirement, such as a finite total number of particles. In 
point of fact, appropriate normalization of Tsallis' distri- 
bution requires a g-parameter greater than 1/3 0]. 

To complete the proof, we now show that Tsallis' equi- 
librium g-distribution || is a natural consequence of the 
-ff 9 -theorem. As happens in the canonical iJ-theorem, 
G q = must be a necessary and sufficient condition for 
equilibrium. Since the integrand appearing in the ex- 
pression of G q cannot be negative, this occur if and only 
if 



ln r /' + hv fx = ln 9* / + U V h 



(16) 



Therefore, the above sum of g-logarithms remains con- 
stant during a collision: it is a summational invariant. 
Only the particles total mass, energy, and momentum 
behave like that ■ Consequently, we must have 



hiq* f = a + d{ ■ v + ci2V ■ v 



(17) 
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where a and 02 are constants and d[ is an arbitrary 
constant vector. By introducing the barycentric velocity, 
u, we may rewrite (|l7|) as 



In,* / = a - j*\v- u\ 



(18) 



with a different set of constants. Taking A q - = e 
and defining 7 = /,J , we obtain a generalized 



•1211/1-8* 



(19) 



JT- 

Maxwell's distribution 

f (v) = A q ,[l-(l-q*) 7 \v- 

where A q * , 7 and u may be functions of the temper- 
ature. The above expression is the general form of the 
q-Maxwellian distribution function |9J . 

Summing up, we have discussed a q-generalization of 
Boltzmann's kinetic equation along the lines of Tsallis 
nonextensive thermostatistics. Our main results followed 
from a slightly modified version of the statistical hypothe- 
ses underlying Boltzmann's approach, incorporating (i) 
the nonextensivity property, explicitly introduced trough 
a new functional form for the local entropy, and (ii) a non- 
factorizable expression for the molecular chaos hypothe- 
sis. Both ingredients were shown to be consistent with 
the standard laws of (microscopic) dynamics. They re- 
duce to the familiar Boltzmann assumptions in the exten- 
sive limit q — ► 1. The usual statistical hypothesis of com- 
pletely uncorrelated colliding molecules seems to be too 
restrictive. It is conceivable that correlations may be rele- 
vant within some scenarios. Here we have provided a sim- 
ple type of correlations that makes sense within Tsallis' 
nonextensive thermostatistics. Other possibilities, also 
leading to Tsallis' distribution (111), are obtained if one 
replaces the function e q {x) in (|9|) by other positive, in- 
creasing function F q (x) such that lim 9 ^i F q (x) = exp(a;). 
Naturally, these (/-generalizations of the molecular chaos 
hypothesis do not settle the profound conceptual issues 
raised by Boltzmann's "Stozssahlansatz" . What we are 
advocating is that Boltzmann's statistical assumptions 
do not encompass all the possibilities allowed by the gen- 
eral principles of mechanics. 

The study of chaotic, low dimensional dissipative dy- 
namical systems has suggested a deep connection be- 
tween Tsallis formalism and multifractals (see and 
references therein). It would be interesting to explore if 
this relationship also holds for Hamiltonian systems of 
large dimensionality and if it does, whether there is any 
connection with the f/g-theorem. 

Finally, we stress that the solutions of the generalized 
Boltzmann equation (]?]) verify the iJ 9 -theorem only if 
q > 0. In that case H q is an increasing function of 
time and the time dependent solutions of (^) evolve ir- 
reversibly towards Tsallis' equilibrium distribution (19). 
These results can be extended to include nonuniform sys- 
tems as well as more general interparticle interactions. 
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